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A  wealth  of  experimental  data  has  shown  that  atmospheric  turbulence  can  be  anisotropic;  in  this  case,  a 
Kolmogorov  spectrum  does  not  describe  well  the  atmospheric  turbulence  statistics.  In  this  paper,  we  show  a  quan¬ 
titative  analysis  of  anisotropic  turbulence  by  using  a  non-Kolmogorov  power  spectrum  with  an  anisotropic  coeffi¬ 
cient.  The  spectrum  we  use  does  not  include  the  inner  and  outer  scales,  it  is  valid  only  inside  the  inertial  subrange, 
and  it  has  a  power-law  slope  that  can  be  different  from  a  Kolmogorov  one.  Using  this  power  spectrum,  in  the  weak 
turbulence  condition,  we  analyze  the  impact  of  the  power-law  variations  a  on  the  long-term  beam  spread  and 
scintillation  index  for  several  anisotropic  coefficient  values  f.  We  consider  only  horizontal  propagation  across 
the  turbulence  cells,  assuming  circular  symmetry  is  maintained  on  the  orthogonal  plane  to  the  propagation  direc¬ 
tion.  We  conclude  that  the  anisotropic  coefficient  influences  both  the  long-term  beam  spread  and  the  scintillation 
index  by  the  factor  fi~a.  ©  2011  Optical  Society  of  America 
OCIS  codes:  010.0010,  010.1300,  010.1330. 


1.  INTRODUCTION 

The  atmosphere’s  statistical  behavior  does  not  always  follow 
the  Kolmogorov  power  spectrum  density  model;  rather  it 
follows  different  power  laws  [1].  Several  experimental  data 
and  measurements  indicate  that  the  turbulence  in  the  upper 
troposphere  and  stratosphere,  or  along  nonhomogeneous 
path,  deviates  from  predictions  of  the  Kolmogorov  model 
[2-6].  Other  experiments  showed  that  atmospheric  turbulence 
in  maritime  environments  can  assume  a  different  statistical 
behavior  with  respect  to  Kolmogorov  [7-13],  In  addition,  an¬ 
isotropy  in  stratospheric  turbulent  inhomogeneities  has  been 
experimentally  investigated  [14-17],  and  laboratory  results 
have  shown  that  turbulence  can  be  anisotropic.  These  studies 
cast  doubt  on  the  correctness  of  the  conventional  assumption 
of  isotropic  turbulence  through  the  entire  atmosphere.  Also, 
researchers  have  collected  data  at  high  altitudes  that  indicate 
the  existence  of  anisotropic  turbulence  and  non-Kolmogorov 
behavior  [18]. 

Kon  [19]  qualitatively  analyzed  the  effect  of  the  anisotropic 
turbulence  on  the  long-term  beam  spread  and  scintillation 
index  using  the  classical  Kolmogorov  power-law  exponent 
a  =  11/3.  In  this  paper,  we  show  a  quantitative  analysis  of  the 
anisotropic  turbulence  by  using  a  non-Kolmogorov  power 
spectrum  with  an  anisotropic  coefficient.  The  spectrum  used 
here  is  based  on  the  spectrum  shown  in  [19];  it  does  not  in¬ 
clude  the  inner  and  outer  scales,  and  it  is  valid  only  inside 
the  inertial  subrange.  Using  this  non-Kolmogorov  spectrum, 
in  weak  turbulence  conditions,  we  analyze  the  impact  of 
exponent  variations  a  on  long-term  beam  spread  and  scintilla¬ 
tion  index  for  several  anisotropic  coefficient  values  c.  We  con¬ 
sider  only  propagation  across  the  turbulence  cells  where 
circular  symmetry  is  maintained  on  the  orthogonal  plane  to 
the  propagation  direction.  In  other  words,  we  suppose  that 
the  anisotropic  turbulence  cells  are  a  solid  generated  by 
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the  rotation  of  an  ellipse  around  its  semiminor  axis,  which 
is  overlapped  with  the  horizontal  propagation  direction.  This 
hypothesis  is  essential  to  maintain  the  symmetry  assumptions 
reported  in  [20]  used  to  calculate  turbulence  statistics  param¬ 
eters,  such  as  scintillation  index.  Although  this  hypothesis  is 
quite  strong,  it  should  be  considered  as  a  first  step  toward  a 
better  understanding  of  a  very  complex  problem  such  as  ani¬ 
sotropic  turbulence.  Finally,  let  us  note  that  often  is  more 
appropriate,  in  the  context  of  concrete  atmospheric  problems, 
to  consider  anisotropic  propagation  along  the  vertical  direc¬ 
tion.  In  such  a  case,  our  results  remain  valid  if  we  take  into 
account  the  change  of  C'/  with  altitude. 

2.  ANISOTROPIC  NON-KOLMOGOROV 
SPECTRUM 

Published  theoretical  papers  [21,22]  were  dedicated  to  the  in¬ 
terpretation  of  experiments  on  the  scattering  of  radio  waves, 
particularly  on  the  nonisotropic  character  of  turbulence  scat¬ 
tered.  Such  papers  used  a  nonisotropic  power  spectrum  in 
agreement  with  experimental  data.  Later,  Kon  [19]  showed  a 
qualitative  study  of  nonisotropic  turbulence  by  using  the  men¬ 
tioned  power  spectrum.  In  this  paper,  we  assume,  as  in  [19], 
that  for  a  laser  beam  propagating  along  a  path  exhibiting  an¬ 
isotropic  turbulence,  the  structure  function  for  the  refractive- 
index  fluctuations  is  given  by 

Dn(r)  =f)-Ci-  (A*2  +  Ag/2  +  A \  (1) 

where  f  is  a  vector  spatial  variable,  y  is  the  power  law,  which 
reduces  to  2/3  in  the  case  of  conventional  Kolmogorov  turbu¬ 
lence.  Here,  ft  is  a  constant  equal  to  unity  when  y  =  2/3,  but 
otherwise  has  units  mrY  1 2/:!  and  c  is  the  anisotropic  coeffi¬ 
cient.  The  corresponding  power-law  spectrum  associated  with 
such  structure  function  takes  the  form  [20] 
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<&„(£,  a)  =  A{a)  ■  Cl  •  e2  ■  [x2  +  c(4  +  4)1  §’  (2) 

where  a  =  y  +  3  is  the  spectral  index  or  power  law,  C„  =  [i  ■ 
C2  is  a  generahzed  structure  parameter  with  units  m~r ,  and 
A  (a)  is  defined  by 

A(a)  =  ^^r(a  -  1)  cosT^V  3  <  a  <  4,  (3) 

and  the  symbol  T(a;)  in  the  last  expression  is  the  ganmia  func¬ 
tion.  When  a  =  11/3  and  t;  =  1,  we  fmd  that  A(  11/3)  =  0.033 
and  the  generalized  power  spectrum  reduces  to  the  conven¬ 
tional  Kolmogorov  spectrum.  Also,  when  the  power  law 
approaches  the  limiting  value  a  =  3,  the  function  A(a)  ap¬ 
proaches  zero.  Consequently,  the  refractive-index  power 
spectral  density  vanishes  in  this  limiting  case. 

3.  LONG-TERM  BEAM  SPREAD 

It  is  well  known  that  the  beam  spot  size  of  a  laser  beam 
propagating  in  turbulence  is  affected  from  two  main  effects: 
beam  spreading  and  beam  wander.  Random  temperature  fluc¬ 
tuations  of  the  atmosphere  yield  random  refractive-index  fluc¬ 
tuations;  therefore,  a  laser  beam  propagating  through  the 
atmosphere  is  randomly  deviated  from  the  direction  of  propa¬ 
gation  (beam  wander),  and  it  is  more  affected  from  beam 
spreading  than  a  diffraction-limited  beam  in  absence  of  turbu¬ 
lence.  Physically,  the  turbulence  acts  like  many  lenses  of  dif¬ 
ferent  size  that  randomly  change  the  effective  optical  path  of 
the  beam.  The  combined  effects  of  beam  wander  and  beam 
spreading  is  called  long-term  beam  spread,  which  represents 
the  effective  beam  spot  size,  and  it  is  used  as  one  of  the  main 
parameters  to  evaluate  the  intensity  profile  along  the  path. 

The  analytical  form  of  long-term  beam  spread  for  a 
Gaussian  beam  is  [20] 

W2e  =  WlT(a)  =  W2-[1  +  T(a)},  (4) 


o-|(a)  =  -8;r2  •  A(a)  ■  -  ■  T ^1  -  ■  sin ■  C\  ■  k3  “  ■  Li. 

(6) 

Note  that,  for  a  =  11/3  and  <;  =  1,  we  obtain  the  particular 
case  of  the  Kolmogorov  spectrum  already  shown  in  [1].  Also, 
for  c  =  1,  Eq.  (5)  reduces  to  isotropic  expression  shown  in  [1]. 

Then,  we  plot  in  Fig.  1  the  long-term  beam  spread  as  a 
function  of  alpha  for  several  c  setting  the  following  param¬ 
eters  (W()  is  the  beam  spot  radius):  L  =  50000  m,  C2  = 
7  ■  10-18  nr“+3,  X  =  1.55 /im,  and  W0  =  0.01  m.  The  value  of  C2 
has  been  chosen  to  simulate  a  high-altitude  path,  where  ani¬ 
sotropic  turbulence  may  be  present. 

We  deduce  from  Fig.  1  that  the  anisotropic  parameter  <; 
influences  the  long-term  beam  spread,  which  is  reduced  if  <; 
assume  high  values.  Also,  we  deduce  from  Fig.  1  that  alpha 
variations  have  less  impact  on  long-term  beam  spread  for  high 
c  values.  These  deductions  can  be  physically  explained  by 
mentioning  the  change  of  curvature  of  the  anisotropic  turbu¬ 
lence  cells  with  respect  to  the  isotropic  case.  Anisotropic  cells 
change  the  focusing  properties  of  the  turbulence;  in  particu¬ 
lar,  a  beam  propagating  along  the  short  axis  of  the  anisotropic 
cells  (z  direction)  will  be  less  deviated  from  the  direction  of 
propagation  because  these  cells  act  as  lenses  with  a  higher 
radius  of  curvature.  Accordingly,  the  long-term  beam  spread 
will  be  reduced,  and,  for  high  values  of  c,  this  anisotropic  ef¬ 
fect  dominates  on  the  power-law  variations. 

We  cannot  deduce  from  this  analysis  how  much  impact  an¬ 
isotropy  has  on  long-term  beam  spread  with  respect  the  inner 
and  outer  scales,  which  are  not  included  in  the  spectrum 
Eq.  (2);  however,  that  could  be  an  interesting  future  task. 

4.  SCINTILLATION  INDEX 

A  laser  beam  propagating  in  turbulence  is  affected  from  an¬ 
other  main  effect:  the  scintillation.  The  scintillation  index  is 
denoted  by 


where  W  is  the  diffraction-limited  spot  size  radius,  and  T(a)  is 
the  term  that  includes  small-scale  beam  spreading  and  beam 
wander  atmospheric  effects  [20]. 

If  we  consider  that  circular  symmetry  is  maintained  on  the 
orthogonal  plane  to  the  propagation  direction,  we  can  use  the 
same  procedure  as  reported  in  [20],  but  using  the  power 
spectrum  Eq.  (2).  We  carry  out 


T(a)  =  4tt2/c2L  •  (J  J  k  ■  4>n{a,  /c)d/cd£ 

-  1  1  K  ■  </>n(a,  k )  exp  ALk^  ^  d<cd/j 


=  0.25  •- 


a 


a-  1 

isotropic  *  ^  ' 


n 

sm  |  a  ■  — 


-i 


■  4(a)- AH- C2 


(5) 


where  f  =  1  - 1,  z  is  the  propagation  distance,  L  is  the  path 
length,  A  =  ^  k  =  i/)n(a,K)  is  the  non-Kolmogorov 

spectrum  Eq.  (2),  and  we  have  defined  a  non-Kolmogorov 
Rytov  variance  as  the  plane  wave  scintillation  index  in  non- 
Kolmogorov  turbulence  [23]: 


(7) 


We  as  a  function  of  alpha  for  several  zeta  values 


Fig.  1.  (Color  online)  Long-term  beam  spread  as  a  function  of  alpha 
for  a  horizontal  path. 
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where  /  is  the  optical  field  irradiance  and  the  brackets  denote 
an  ensemble,  or  long  spatial  distance,  average. 

In  our  analysis,  as  reported  in  [22],  we  include  aperture 
averaging  effects  of  the  receiver  aperture  of  diameter  DG, 
so  we  carry  out  the  flux  variance  in  the  plane  of  the  detector 
at  a  short  distance  Lf  behind  the  collecting  lens.  We  illustrate 
such  a  system  in  Fig.  2. 


Collecting  lens 


Fig.  2.  Propagation  geometry  for  a  laser  beam  originating  at  distance 
L  to  the  left  of  a  thin  Gaussian  lens. 


A.  PLANE  WAVE  MODEL 

The  procedure  adopted  in  this  paper  is  the  same  as  discussed 
in  [20]  for  the  standard  Kolmogorov  spectrum,  but  here  we 
use  a  non-Kolmogorov  anisotropic  spectrum  Eq.  (2).  Our  ana¬ 
lysis  for  the  plane  wave  model  leads  to 


where  2F'1(a,  6;  c;«)  is  the  hypergeometric  function  that  is 
given  by 


2  F^a^c^z)  =  J2 


(a)n-(b)n  zn 
(c)n  n\ ' 


M  <  °°>  (9) 
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Note  that  if  we  consider  a  point  receiver  ( DG  =  0)  no  aperture 
averaging  is  present  and  o’|piane(«,flG)  =  Sf(a),  which  is  the 
non-Kolmogorov  Rytov  variance  [Eq.  (6)].  Also,  for  c  =  1, 
Eq.  (8a)  reduces  to  the  isotropic  expression  shown  in  [1]. 

We  plot  in  Fig.  3  c-ypl,me(a,  I)G)  as  a  function  of  alpha  for 
several  c  for  a  particular  high-altitude  horizontal  path.  We 
use  the  following  set  of  parameters: 

L  =  50000 m;  C\  ~  7  ■  10“ 18  nr“+3; 

/l  =  1.55/un;  OG  =  0.1m. 


B.  SPHERICAL  WAVE  MODEL 

Similar  to  the  plane  wave  model  case,  the  scintillation  analysis 
for  the  spherical  wave  model  leads  to 


where 


(«)» 


T(a  +  n) 


(10) 


We  plot  in  Fig.  4  fffSphericai(a.-DG)  as  a  function  of  alpha  for 
several  c  for  a  particular  high-altitude  horizontal  path,  taking 
the  same  set  of  parameters  as  the  plane  wave  case. 


C.  GAUSSIAN  BEAM  WAVE  MODEL 

To  describe  the  beam  characteristics  at  the  input  plane  and 
the  front  plane  of  the  Gaussian  lens,  we  use  two  sets  of  non- 
dimensional  beam  parameters.  We  assume  the  transmitted 
Gaussian  beam  at  the  input  plane  has  finite  radius  W0  and 
phase  front  radius  of  curvature  given  by  F0.  Thus,  we  have 


J I  spherical 


(a.DG)  =  8k2  ■  k2  ■  L  ■  Rei^J  J  k  ■  (k,  a)  ■  exp  f-  ^ — — 
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Cn2=7e-18,L=50000m,Dg=0.1m 


Fig.  3.  (Color  online)  Scintillation  index  for  the  plane  wave  as  a 
function  of  alpha  for  several  anisotropic  coefficient  values. 


x  io'3  Spherical  wave  setting  Cn2=7e-18,L=50000m,Dg=0.1  m 


Fig.  4.  (Color  online)  Scintillation  index  for  the  spherical  wave  as  a 
function  of  alpha  for  several  anisotropic  coefficient  values. 


L  2 L 

at  the  transmitter  (z  =  0):  0O  =  1  -  — ,  A0  = - ^  >  (11) 

F  o  kWf. 


at  the  Gaussian  lens  (z  =  L):  0t  =  ,^° 


0g  +  Ag' 


©g+Ag’ 

(12) 


Following  the  same  procedure  discussed  in  [1]  for  the 
standard  Kolmogorov  spectrum,  but  here  using  the  non- 
Kolmogorov  anisotropic  spectrum  Eq.  (2),  our  analysis  for 
a  Gaussian  beam  leads  to  (neglecting  beam  wander  effects) 


of  the  anisotropy  coefficient  g  decrease  the  scintillation  index 
and,  at  the  same  time,  reduce  the  impact  of  a  variation.  These 
deductions  can  be  physically  explained  by  mentioning  the 
change  of  curvature  of  the  anisotropic  turbulence  cells  with 
respect  to  the  isotropic  case.  Turbulence  cells  act  as  lenses 
with  a  higher  radius  of  curvature,  leading  to  a  reduction  of 
amplitude  fluctuations  and,  consequently,  scintillation  index. 
As  for  long-term  beam  spread,  we  cannot  deduce  from  this 
analysis  how  much  impact  anisotropy  has  on  the  scintillation 
index  with  respect  to  the  inner  and  outer  scales.  For  g  =  1, 
Eq.  (13)  reduces  to  the  isotropic  expression  shown  in  [24]. 
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We  plot  in  Fig.  5  ^Gaussian  F(;)  as  a  function  of  alpha  for 
several  g  for  a  collimated  beam  (0O  =  1),  taking  the  same 
set  of  parameters  previously  used. 

We  deduce  from  Figs.  3-5  that  the  anisotropic  coefficient 
influences  the  scintillation  index.  In  particular,  high  values 


5.  ANISOTROPIC  FACTOR 

The  main  result  deduced  from  all  the  previously  reported  ex¬ 
pressions  is  the  multiplicative  factor  c2~",  which  essentially 
introduces  a  turbulence  rescaling  due  to  anisotropic  and 
non-Kolmogorov  turbulence. 

We  plot  the  anisotropic  factor  c2~"  as  a  function  of  a  for 
several  g  values.  Results  are  shown  in  Fig.  6. 


Cn2=76v18,L=50000m,Wo=0.01m,Dg=0.1m 


Fig.  5.  (Color  online)  Scintillation  index  for  the  Gaussian  beam  as  a 
function  of  alpha  for  several  anisotropic  coefficient  values. 
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Fig.  6.  (Color  online)  Anisotropic  factor  as  a  function  of  alpha  for 
several  zeta  values. 


We  deduce  from  Fig.  6  that  the  anisotropic  factor  decreases 
as  the  alpha  value  increases  for  every  zeta  value.  Also,  high 
zeta  values  reduce  the  anisotropic  factor  and,  consequently, 
the  atmospheric  turbulence  parameters  (scintillation  index 
and  long-term  beam  spread). 

In  this  paper,  we  supposed  that  the  anisotropic  factor  has 
the  same  effect  on  all  the  turbulence  scales.  However,  in 
stable  atmospheric  boundary  layers,  isotropy  probably  pre¬ 
vails  at  small  scales.  Therefore,  in  this  case,  we  should  distin¬ 
guish  small  scales  from  other  scales  by  using  a  more  complete 
power  spectrum  with  a  filter  function.  Finally,  we  cannot  de¬ 
duce  from  Fig.  6  how  much  the  anisotropic  factor  is  significant 
in  practice,  given  all  the  inaccuracies  in  estimating  structure 
function  parameters,  outer  scales,  etc.,  that  occur.  Future  ex¬ 
periments  will  give  us  an  idea  about  the  validity  of  this  theo¬ 
retical  model. 

6.  DISCUSSION 

In  this  paper,  we  used  a  non-Kolmogorov  refractive-index 
structure  function  and  a  non-Kolmogorov  power  spectrum 
with  an  anisotropic  coefficient  q  to  describe  anisotropic  tur¬ 
bulence.  We  assumed  that  the  circular  symmetry  is  main¬ 
tained  on  the  orthogonal  plane  to  the  propagation  direction 
z,  and  we  analyzed,  in  weak  turbulence  condition,  long-term 
beam  spread  and  scintillation  index  for  a  horizontal  path  and 
for  several  anisotropic  coefficient  values  q.  We  concluded  that 
the  anisotropic  parameter  influences  both  the  long-tenn  beam 
spread  and  the  scintillation  index  by  the  factor  <^~a.  Therefore, 
a  high  value  of  q  decreases  both  the  long-term  beam  spread 
and  the  scintillation  index.  Finally,  we  showed  that  power- 
law  a  variations  have  less  impact  on  both  the  long-term  beam 
spread  and  the  scintillation  index  as  q  increase. 

These  conclusions  can  be  physically  explained  by  mention¬ 
ing  the  change  of  curvature  of  the  anisotropic  turbulence  cells 
with  respect  to  the  isotropic  case.  Anisotropic  cells  change 
the  focusing  properties  of  the  turbulence;  in  particular,  a  laser 
beam  propagating  along  the  short  axis  of  the  cells  ( z  direc¬ 
tion)  will  be  less  deviated  from  the  direction  of  propagation 
because  these  cells  act  as  lenses  with  a  higher  radius  of 
curvature.  Accordingly,  the  long-term  beam  spread  will  be  re¬ 
duced,  and,  for  high  values  of  q,  this  anisotropic  effect  dom¬ 


inates  on  the  a  variations.  For  the  same  reason,  amplitude 
fluctuations  decrease  as  compared  with  the  isotropic  case, 
leading  to  a  reduction  of  the  scintillation  index. 

Unfortunately,  we  cannot  deduce  from  this  analysis  how 
much  impact  anisotropy  has  on  both  the  long-term  beam 
spread  and  the  scintillation  index  with  respect  the  inner  and 
outer  scales.  Also,  we  supposed  that  the  anisotropic  factor  has 
the  same  effect  on  all  the  turbulence  scales;  however,  in  stable 
atmospheric  boundary  layers,  isotropy  probably  prevails  at 
small  scales.  Therefore,  in  this  case,  we  should  distinguish 
small  scales  from  other  scales  by  using  a  more  complete 
power  spectrum,  which  should  include  also  a  filter  function. 
Finally,  we  are  aware  that  this  anisotropic  turbulence  spec¬ 
trum  needs  to  be  validated  by  comparing  theoretical  results 
with  experimental  data.  Therefore,  this  paper  should  be  con¬ 
sidered  the  first  step  of  a  more  complete  analysis  that  should 
be  based  on  new  experiments. 
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